Synchronization interfaces and overlapping communities in complex networks 
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We show that a complex network of phase oscillators may display interfaces between domains 
(clusters) of synchronized oscillations. The emergence and dynamics of these interfaces are studied 
in the general framework of interacting phase oscillators composed of either dynamical domains 
(influenced by different forcing processes), or structural domains (modular networks). The obtained 
results allow to give a functional definition of overlapping structures in modular networks, and 
suggest a practical method to identify them. As a result, our algorithm could detect information on 
both single overlapping nodes and overlapping clusters. 



The functioning of many natural (biological, neural, 
chemical) or artificial (technological) networks displays 
coordination of parallel tasks [l| . This phenomenon may 
be represented as the interplay between two simultane- 
ous processes. The first (involving most of the network 
nodes) leads to the emergence of organized clusters (or 
moduli, or cohesive subgroups), where nodes in the same 
cluster adjust their dynamics into a common (synchro- 
nized) behavior to enhance the performance of a specific 
task. The second process (involving just few nodes of 
the graph) is to form interfaces (or overlapping struc- 
tures) between the moduli, that are responsible for the 
coordination between the different tasks. 

A practical example is the brain functioning, wherein 
the vision of the images from the right and left eye can be 
represented as the emergence of a collective behavior in 
different, close-by areas of neurons, whereas perception 
(as e.g. feature binding) implies the coordination of the 
two emergent dynamics |!2:|]. Another example in social 
science is, for instance, consensus formation in polarized 
elections (as e.g. between a left-wing and a right-wing 
candidate for the presidency) , where usually most of cit- 
izens cast their votes following a consistent opinion (they 
always vote for the left-wing or the right-wing party), 
whereas usually a small fraction (mostly being respon- 
sible for the final election's outcome) alternate in time 
their preference depending on their actual opinion [3|. 

In this Letter, we present the first evidence that, un- 
der the presence of different functional (synchronized) 
clusters, interfaces appear and show a novel specific dy- 
namical behavior. This finding enables us to develop an 
algorithm for extracting information on the overlapping 
structure in a modular network. Indeed, the study of 
separate modular structures J4| and synchronization [5] 
in complex graphs has not so far unraveled the crucial 
point concerning the role of synchronization interfaces 
and its usefulness in detecting overlapping communities. 

In the following, without lack of generality, we study 
networks consisting of two domains of interacting phase 
oscillators (each one synchronously evolving at a differ- 
ent frequency) , where the nature of the two different fre- 
quency domains is the result of either a dynamical pro- 



cess (influenced by different forcing processes) or a struc- 
tural design (modular network). Under these conditions, 
at each time, most of the oscillators will contribute to 
the synchronous behavior of the two clusters, whereas a 
few nodes will find themselves in a frustrating situation 
of having to decide how to behave as a consequence of 
the contrasting inputs received by the two clusters. 

Let us first describe the case of a synchronization in- 
terface emerging from a dynamical process in a generic 
random graph G of A^ coupled oscillators, whose origi- 
nal frequencies {loi} are randomly drawn from a uniform 
distribution in the interval 0.5 ±0.25, subjected simulta- 
neously to an internal bidirectional coupling and an ex- 
ternal pace-making unidirectional forcing. The network 
dynamics is described by: 






(1) 



yf^i ~r f^p 



where dots denote temporal derivatives, ki is the degree 
of the i*'* oscillator, 0p. is the instantaneous phase of a 
forcing oscillator having fcp. unidirectional connections, 
d and dp are coupling strengths, and {fly} are either 1 
or depending on whether or not a link exists between 
node i and node j . 

In our simulation, we study a network G which con- 
sists of A^ = 200 phase oscillators arranged in an Erdos- 
Renyi configuration [6|. Initially, we set dp — Q and 
kp. — 0\/i, and we choose d = 0.1 such that G features an 
unsynchronized motion. Next, we arbitrarily divide the 
nodes into two groups: nodes from i = 1, • • • , 100 (from 
i = 101, • • • , 200) are assigned to the community A (to 
the community B). We introduce two pacemakers of fre- 
quencies ujp^ and ujpj^ , and connect the nodes in the first 
group (in the second group) with the first (second) pace- 
maker. This implies in Eq. ([1]) to set </>p. = <^„^ = w„^t 
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Opg t) for all the nodes in A {B) . In order to 
assign to each node its kp. links with the pacemaker, we 
start at io = the evolution of Eq. ^ from random ini- 
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FIG. 1: (Color online) (a-c) Instantaneous frequencies (j}i{t) 
of each one of the 100 oscillators in A (red lines) [in B, (blue 
lines)] vs. time, obtained from simulation of Eq. ^ for (a) 
d = 3.25, (b) d = 5.75, and (c) d = 9.75. (d) Log-log plot 
of the period To of the switching process in the interface vs. 
uja, for d — 3.5. The solid line represents a linear fit with 
To ~ 120/a;A- Each point is the average of 5 independent 
realizations. Other parameters reported in the text. 



tial conditions in the unforced case {kp. — Vi), and add 
at later times ti = to + lAt links between nodes of the two 
communities and the pacemakers. Precisely, at each time 
ti , the pace-maker (pp^ {4>pb ) forms a connection with 
that node j in A (B) whose instantaneous phase at time 
ti corresponds to the minimum miuj \S — A9j mod 2tt\ , 
with A9j = 4'ji'ti) ~ <t>p{ti) 1 ^ ^ (0, 27r) a suitable param- 
eter, and <pp{ti) = <i>uA.ii) (0p(*/) = 'Ppsiii)) for those 
nodes in A (B). Ref. [3] demonstrated that, by operating 
this attachment over a given time interval, the resulting 
dynamics of any arbitrary network of oscillators can be 
entrained to any arbitrary frequency ujp. 



By selecting Wp^ = 0.7, ujpg — 0.3, J2i=i ^Pi — 2,000 
links to the pacemakers (1,000 to entrain the nodes in A 
and the other 1,000 to entrain those in B), and dp = 1, 
the dynamics displays two large communities of entrained 
oscillators, densely intermingled by the connectivity of G. 
The situation is depicted in Fig.[T]Ja) where we report the 
instantaneous frequency of each oscillator in G (averaged 
over a small window to smooth fluctuations) as a func- 
tion of time. We observe that most of the nodes in A (in 
B) have a constant frequency (that of the corresponding 
pacemaker) , whereas the few nodes belonging to the syn- 
chronization interface feature a frequency which is oscil- 
lating around the mean value of that of the two communi- 
ties u! — {ujp^ +ujpg)/2. Furthermore, Fig. [TJd) gives ev- 
idence that the period Tq of the switching process in the 
interface is inversely proportional to wa = (i-^pA ~'^pb )/2- 

Notice that this switching mechanism is the result of 
the competition of two conflicting processes: the syn- 
chronization within G (controlled by the parameter d) 



that would lead the whole network to exhibit a unique 
frequency, and the forcing of the two pacemakers (con- 
trolled by dp) that would try to separate the nodes into 
two clusters of entrained oscillators. In order to quantify 
this competition and to describe the size of the synchro- 
nization interface we fix dp = 1 and gradually increase d. 
The results are shown in Figs. [Ija-c) with (a) d = 3.25, 
(b) d = 5.75, and (c) d = 9.75. For intermediate coupling 
Fig. [ijb), it is observed that the interface grows in size, 
recruiting more and more nodes out of the two commu- 
nities into the oscillating mode. Due to the presence of 
the two forcing pacemakers, this interface is organized in 
an oscillating mode rather than in a constant frequency 
mode, as it would occur in the classical Kuramoto regime 
Q. As the coupling d is further increased, almost the en- 
tire system of oscillators eventually participates into this 
interface oscillating mode, as seen in Fig. [Ijc). Notice 
that similar switching dynamics was observed before in 
the case of a chain of oscillators subjected to two forcing 
frequencies applied to the two ends of the chain @ . 

In order to give an analytical insight to our findings, 
let us consider the simple case of three interacting phase 
oscillators described by: (pi — uji + Ki sin((/)3 — 0i) ]4>2 — 
W2+^2sin((^3-(/)2) ;03 = W3-|-ii^[sin(0l-03)-|-sin(02- 
(/l3)]. Here, 03 is the phase of an oscillator (with natu- 
ral frequency W3) receiving a simultaneous coupling from 
two other oscillators at natural frequencies wi 7^ W2 7^ W3, 
and Ki, K2 <C K are coupling constants. The oscillators 
1 and 2 model the two functional (frequency) domains A 
and B, where the nature of the two different frequency 
domains is from either a dynamical process or a struc- 
tural design. These two domains of synchronous oscil- 
lators simultaneously interact with the small group of 
nodes in the interface (modeled by the third oscillator), 
so that we can reasonably assume Ki = K2 = 0, and 
consequently 0i_2 — 0Ji,2t- 

The resulting equation for 03 



03 = CJ3 + 2A sm ( ;:; 1 - 03 ) cos ( 1 (2) 



yields an analytic solution 
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(3) 



for 0^3 = Lu, and where ^3 = 03 — cut, u ~ (wi + uj2)/2, 
LUj\ — (u}\ — L02)/1, and A,B suitable parameters. In 
good agreement with the results of Fig. [TJd), Eq. ([3]) 
predicts that the instantaneous frequency 03 (once ^3 is 
selected to be the mean frequency of the two forcing clus- 
ters) will oscillate around Q with a period Tq inversely 
proportional to the difference in the frequencies of the 
two forcing clusters. Notice that, as K increases (i.e. in 
the strong coupling regime), 6^ will progressively van- 
ish, which means that 03 = a) or, in other words, the 
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FIG. 2: (a) Instantaneous frequencies (j}i{t) vs. time from 
simulation of Eq. ^ with d = 0.1 (other parameters and 
stipulations are reported in the text). Squares, diamonds, 
and full circles represent respectively nodes belonging to A, 
B and O. (b) Log-log plot of the switching period To of 
the oscillations in the frequency of the nodes in O vs. the 
frequency difference u^. The solid line represents a linear fit 
with slope -1.002 ± 0.0069. 



frequency of the oscillator will be locked to the mean 
frequency of the two giant clusters. 

While so far we have considered the behavior of in- 
terfaces as the result of the competition of dynamical 
domains, we now move to describe the case of synchro- 
nization interfaces under the competition of structural 
domains in modular graphs in the absence of the forc- 
ing [dp — and fcp. = Vi in Eq. ([T])]. For this pur- 
pose, we construct the adjacency matrix of G by con- 
sidering two large communities {A, B), each one formed 
by 50 densely and randomly connected nodes (the av- 
erage degree in the same community is 16), which are 
overlapped by a small community O made of a complete 
graph of 5 nodes(3 links to A {B)) that form symmetric 
connections to nodes in both A and B. Each of the 105 
nodes of G is associated with a phase oscillator obeying 
Eq. ([T]), which is integrated with a Runge-Kutta-Fehlberg 
method, for an initial distribution of frequencies such 
that nodes in A (B) [i.e. nodes from i = 1 {i — 51) 
to i = 50 (i = 100)] have natural frequencies uniformly 
distributed in the range 0.25 ± 0.25 (0.5 ± 0.25), while 
nodes in O have frequencies uniformly distributed in the 
range 0.375 ± 0.05) (i.e. around the mean frequency of 
the two distributions). 

Figure [2] shows that all oscillators in the communi- 
ties A and B behave synchronously with an almost con- 
stant frequency in time (well approximating the mean 
of the original frequency distribution), while all oscilla- 
tors in O constitute the synchronization interface and, 
as so, display an instantaneous frequency oscillating in 
time around the mean value of the two frequencies in 
the two clusters [Fig. [D[a)]. Similarly, Fig. ^h) shows 
the good agreement between the analytical prediction of 
Eq. ([3]) and the numerical results concerning the scaling 
of the period of the frequency oscillations of O with the 
frequency difference between the two communities A and 
B. 

It is important to remark that all other analytical pre- 
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FIG. 3: Effects of coupling strength (a-b) and asymmetry (c). 
(a) d = 0.5; (b) d = 0.95. (c) d = 0.1, but 2 (5) links of the 
nodes in the interface are to nodes of cluster A (B). As in 
Fig. [5] squares, diamonds and full circles refer respectively to 
nodes belonging to A, B and O. 



dictions of Eq. ^ concerning the behavior of the syn- 
chronization interface in the large coupling regime can be 
confirmed in our simulations. Indeed, looking at Fig.[3][a- 
b) one realizes that the effect of increasing the coupling 
strength d is to lock almost always the synchronized fre- 
quency of the interface to the mean of the frequencies 
of the two communities, with persistent events of short 
shootings to even larger (or lower) frequencies. 

What described so far refers to the specific case in 
which initially the nodes of the interface were prepared to 
have the same number of links to nodes of the two com- 
munities. It is therefore interesting to ask what happens 
when the nodes in O are asymmetrically connected to A 
and B. While the low coupling regime does not substan- 
tially differ from the symmetrical case, the high coupling 
regime [illustrated in Fig. [Sj^c)] features frequency oscil- 
lations of the nodes in O that are biased toward that 
community to which the nodes have more connections. 
In fact, it can be proved that if the nodes in O have kf 
{kf) connections to nodes in A (B), they tend to lock to 



kfiUA+kfi^E 



k^+k^ 



where 



the weighted mean frequency uj^, 
^A (wb) is the average frequency of the nodes in A (B). 
The overall scenario reported above suggests a practi- 
cal way to detect overlapping communities (Fig.JD^a)) in 
generic modular networks. It is essential to remark that 
most of the definitions of network communities proposed 
so far led essentially to a graph partition into compo- 
nents, such that a given node belongs to and only to one 
of the components of the partition [4]. The possibility, 
instead, that two components of a partition can have an 
overlapping set of nodes has been recently investigated 
by means of topological arguments [lO] . The novelty of 
our approach is in introducing a functional concept of 
overlapping structures, that are defined in relationship to 
dynamical response of the network as a whole. Namely, 
as far as synchronized behavior of phase oscillators is con- 
cerned, we define an overlapping structure as the set of 
nodes which, instead of following the constant frequency 
of one of the two domains, balance their instantaneous 
frequencies in between these two, and as so, they cannot 
be considered as a functional part of any single domain. 
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FIG. 4: (a) Illustration of the concept of overlapping struc- 
ture, where the overlapping region of two circles (commu- 
nities) represents the overlapping community, (b) C (see 
text for definition) vs. node index in G. Overlapping nodes 
1,2,3(101,102,103), labeled respectively with triangle, dia- 
mond and star, have [1, 5], [2, 4], [3, 3] connections with nodes 
of \B, A]{[A, B]), while 4 - 100(104 - 200) labeled with 
squares belong functionally to A (B) and d = 0.15 (c) Over- 
lapping clusters each consisting of complete graphs of nodes 
1 — 4(101 — 104), labeled with stars, have symmetrical connec- 
tions with A and B, while nodes 5 - 100(105 - 200), labeled 
with squares, belong functionally to A (B), here d = 0.1. 



As we will see below, this definition allows detection of 
further information including single overlapping nodes, 
as compared to previous studies. 

To illustrate this idea, we construct a network G made 
of two large moduli {A and B of 100 nodes each), where 
the majority of nodes forms random connections (with 
average degree 15) with only elements of the same com- 
munity, while only very few nodes form links with nodes 
of both communities. Precisely, we call kf (kf) the to- 
tal number of links these few nodes form with nodes in 
A (B), and define a quantity to evaluate the degree of 
overlapping of these nodes as the ratio Di = kf/kf. 
When the node is perfectly overlapped between two clus- 
ters, Di — 1. Under these conditions, Eq. ^ is simu- 
lated for dp ~ 0, kp. — Vi, the set {tOi} drawn from 
a Gaussian distribution with standard deviation 0.1 and 
mean value 0.3 (0.6) for nodes belonging to community 
A {B). To identify those overlapping nodes, the quantity 
Ci = sgn[(l)i(t) — uj] niint{\(f>i(t) — lu\} (w is the mean of 
the two averaged frequencies assigned to the two commu- 



nities) is introduced, which allows to monitor how close 
in time the dynamics of a node gets to oj. When the 
node is in the dynamics of fully synchronization inter- 
face, Ci — 0. Therefore, the closer Di gets to 1, the 
closer Ci is expected to get to 0. The results are shown 
in Fig. [3l^b,c) for two different arrangements of the over- 
lapping community: overlapping nodes [Fig. \Mjo)] and 
overlapping clusters [Fig. |D^c)] which have symmetrical 
connections to two clusters. In both cases, two large syn- 
chronized clusters are identified very far from the over- 
lapping synchronization, corresponding to those nodes 
performing distinct tasks and already classified by the 
structural partition. At the same time, the dynamical 
evolution manifests a group of nodes whose dynamics is 
located significantly out of the two main clusters (thus 
identifying the overlapping community). In Fig. (H^b), 
each overlapping node gives rise to a different value of 
Ci in correspondence to its specific degree of overlapping 
Di. Namely, the node with Di — 1 yields Ci — 0. On 
the contrary, in Fig. \M,c), all the nodes inside the over- 
lapping cluster are identified as a whole and feature the 
same value of Ci = due to symmetrical connection. 

In conclusion, we have analytically and numerically de- 
scribed the behavior of synchronization interfaces in both 
random and modular complex networks of phase oscilla- 
tors. Wc identify the overlapping structure (nodes or 
clusters) in the modular network, using for the first time 
a functional approach. Therefore, our study is of prac- 
tical relevance for applications to large biological, neu- 
ral, chemical, social and technological modular networks, 
where, besides detection of the structure of overlapping 
communities, one is interested in inspecting the coordi- 
nating role of such overlapping communities in the func- 
tional performance of the graph. 
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